TU-463 
February 1, 2008 



X 



Gravitino Decay into a Neutrino and a Sneutrino in 

the Inflationary Universe 

M. Kawasaki 

Institute for Cosmic Ray Research, The University of Tokyo, Tanashi 188, Japan 



^3 . and 
< 



T. Moroi fl 

Department of Physics, Tohoku University, Sendai 980, Japan 



(N 

cn 

00 

o 

■ Abstract 

I Gravitino produced in the inflationary universe are studied. When the gravitino 

decays into a neutrino and a sneutrino, the emitted high energy neutrinos scatter off the 
background neutrinos and produce charged leptons (mainly electrons and positrons), 
(— I I which cause the electro-magnetic cascades and produce many soft photons. We obtain 

the spectra of the high energy neutrinos as well as the spectrum of the high energy 
photon by integrating a set of Boltzmann equations. Requiring these photons should 
not alter the abundances of the light elements (D, ^He, ^He) in the universe, we can 
I set the stringent upperbound on the reheating temperature after the inflation. We flnd 

that Tr < (10^° - 10^2^GeV for my2 ~ (lOOGeV - ITeV), which is more stringent 
than the constraints in the previous works. 



-'^ Fellow of the Japan Society for the Promotion of Science. 



1 Introduction 



In models based on supergravity [|^, it has been pointed out that there exist some 
hght particles whose interactions are suppressed by powers of the gravitational scale M = 
Mpi/^/Sn ~ 2.4 X lO^^GeV. Such particles have nothing to do with collider experiments, but 
may affect the standard scenario of big-bang cosmology ^, ^. The gravitino, which is the 
gauge field associated with local supersymmetry (SUSY), is one of the weakly interacting 
particles in supergravity models, and we expect that the mass of the gravitino is of order of 
typical SUSY-breaking scale. 

Since the gravitino has only gravitational interaction with other particles, its lifetime is 
very long even if it is unstable. In standard cosmology, the gravitino with mass 7713/2 ~ 
O(100GeV - lOTeV) is excluded since it decays after the big-bang nucleosynthesis (BBN) 
and produces an unacceptable amount of entropy, which conflicts with the predictions of 
BBN [|] . We may avoid this constraint by assuming the inflationary cosmology, in which the 
initial abundance of the gravitino is diluted by the exponential expansion of the universe |0. 
However gravitinos are reproduced by scattering processes off the thermal radiation after 
the universe has been reheated. Since the number density of the secondary gravitino is 
proportional to the reheating temperature, upperbound on the reheating temperature should 
be imposed not to overproduce gravitinos - [|l^. 

In the recent paper |]16|, we consider the unstable gravitino which decays into photon (7) 
and photino (7) and derive upperbound on the reheating temperature Tr (Tr < lO^GeV 
for ~ lOOGeV, Tr < lO^GeV for ~ ITeV) by requiring the produced high energy 
photons should not alter the abundances of light elements synthesized in BBN. The above 
constraint seems to be very stringent since such low reheating temperature requires very 
small decay rate of the infiaton field. For example, in chaotic inflation with a infiaton whose 
interactions are suppressed by gravitational scale M, the infiaton decay rate is expected to 
be Tinf ~ m^^jr/Mpi with niinf being the infiaton mass, and hence the reheating temperature 
is estimated asilTT 



Tr ~ O.l^TinfMpi ~ lO^GeV, (1) 

requiring that the infiaton field should produce the density perturbations observed by COBE 
K„;~10i3GeV) 0. 

However the constraints might become much weaker when the gravitino decays only into 
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weakly interacting particles. In the particle content of the minimal SUSY standard model, 
the only candidates are neutrino and sneutrino. 

In this letter, we assume that the gravitino decays only into a neutrino and a sneu- 
trino and derive constraints on the reheating temperature. In this case, the emitted high 
energy neutrino may scatter off the background neutrino and produce an electron-positron 
(or muon-anti-muon) pair, which then produces many soft photons through electro-magnetic 
cascade processes and destruct the light elements. Since the interaction between high energy 
neutrinos and background neutrinos is weak, it seems that the destruction of the light ele- 



ments is not efficient. In fact, Gratsias, Scherrer and Spergel showed that the constraint 
is not so stringent for the case that the gravitino decays into a neutrino and a sneutrino. 
However the previous analysis seems to be incomplete in a couple of points. First Grat- 



sias et.al. totally neglect the secondary high energy neutrinos which are produced by 
neutrino-neutrino scattering. The effect of the secondary neutrino may be important for 
heavy gravitino (^3/2 > ITeV). Second, they only study the case where the destruction 
of ^He results in the overproduction of ('^He -|- D). However, for the heavy gravitino which 
decays in early stage of the BBN, the destruction of D is more important since the electro- 
magnetic cascade process is so efficient that the energy of soft photons becomes less than 
the threshold of ^He. Furthermore, as pointed out in ref. the previous estimation of the 
gravitino production in the reheating epoch is underestimated. Those effects which are not 



taken into account in ref.|l3| may lead to more stringent constraint on the reheating tem- 
perature. Therefore, in this letter, we reexamine the decay of the gravitino into a neutrino 
and a sneutrino (^/^^ —>■ u + u) with taking all relevant effect into account, and as a result, 
we have obtained more stringent upperbound on the reheating temperature. 

2 Neutrino-Neutrino Interaction 

The high energy neutrinos (z/) produced in the gravitino decay scatter off the thermal 
neutrino (ub) in the background by the following processes; 

T^i + T^i,b T^i + l'i, (2) 

+ l'i,b 1^1 + T^i, (3) 

i^i + i^i,b + (4) 

+ yj^h Uj, (5) 



^i + ^3,b J^i + T^j, (6) 

J^i + e~" + e+, (7) 

+ T>i,h H~ + 1^'^: (8) 

where index i and j represent e, /i and r with i j. All the amplitude square \A4\^ in these 
reactions take the form given by 

= 32G| {a(pp')' + KPQf + <P<l'f + dm\pp')} , (9) 

where Gj? ~ 1.17 x 10~^GeV~^ is the Fermi constant, the coefficients a - d depend on the 
individual reaction, p and p' are the initial momenta of high energy neutrino and background 
neutrino, q and q' are the final momenta, and m represents the mass of the fermion the in 
final state. Then one can obtain the Boltzmann equations describing the evolution of the 
spectra for the high energy neutrinos; 
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dt Sn 



POO _ _ _ _ 

X / dE,Ej,{E,) 
Jo 

^ ^ / 1 1 \ 



\ / i/i+i'i-*e +e+ \ / i/i+i/i^n 

+ ^^n3/2S{El - m3/2/2) 

6r3/2 

+KH^^^^ - 2HUEl)El, (10) 
where H is the expansion rate of the universe, coefficients a - c are given by 

(lout = 4, bout = 0, Cout = 13, (11) 

(^in,ii 6, biji^ii 9, Cin,ii H) (1^) 

and the spectrum of the background neutrino fi, can be written as 
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with being the neutrino temperature. The last two terms of right-hand side of eq.(|T0D 
represent the effect of the cosmic expansion. The formula for charged lepton pair creation 
is obtained as 



dt 
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Coefficients for the charged lepton production processes are given by 



a 



'-^Ui+Ui^L +1+ 



d 



■ui+Pi^r +1+ 



a 



d 



•vi+Pi^lJ+lj 



= 0, 

= IQsm^Owi 

= 16sin'^6'^y, 

= IQsm^Ow^ 

= 0, 

= IQsin^Owi 

= (16 sin^ 9w-2 

= 16 sin"^ 6^14/ — 8 sin^ 



(15) 



(16) 



(17) 



(18) 
(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 



with 6w being the Weinberg angle and i ^ j- The gravitino decay rate r3/2 and the number 
density 113/2 of the gravitino are given by 



"^3/2 

1927rM 



1 - 



(26) 



anc 



0.0232 log 



R 



g-i3/2t (27) 



.lOiOQeV. 

where m^;. is the mass of the sneutrino of z-th generation and is the photon number 
density. In the following analysis, we assume all the sneutrino masses are degenerate for 
simplicity. 

In Fig.|l] the spectra of three types of neutrinos at T = leV are shown for 171^/2 = 
lOOGeV, ITeV and lOTeV. As one can see, the differences among three types of neutrinos 
are very small since the differences only come from the charged lepton pair creation processes 
which are sub-dominant compared with the neutrino-neutrino scattering processes. For the 
heavy gravitino {1713/2 > ITeV) the decay occurs early in the universe when the neutrino 
temperature is so high that the u ~ ui, scatterings is effective. (Notice that the cross section 
is roughly given by (av) ~ G'^pEyTy.) Therefore high energy neutrinos emitted from the 
gravitinos scatter off the background neutrinos frequently and produce many secondary 
neutrinos. On the other hand, the emitted neutrinos seldom interact with the background 
neutrinos if the gravitino mass is small (^3/2 < lOOGeV) and the effect of the neutrino- 
neutrino scattering is negligible. 

3 High Energy Photon Spectrum 

The primary and secondary high energy neutrinos scatter off the background neutrinos 
and produce charged lepton pairs. Then the charged leptons induce the electro-magnetic 
cascade processes. In the same procedure as ref. we obtain the high energy photon 



and electron spectra, taking the following radiative processes into account. (I) High energy 

photon with energy can become electron-positron pair by scattering off background photon 

if the energy of the background photon is larger than m\/e^ where irie is the electron mass. 

We call this process double photon pair creation. For sufficiently high energy photons, this 

is the dominant process since the cross section or the number density of target is much 

larger than other processes. Numerical calculation shows that this process determines the 

shape of the spectrum of high energy photon for > m1/22T . (II) Below the effective 

threshold of double photon pair creation, high energy photons lose their energy by photon- 

^In the recent article, Fischler p9[ | have proposed a new mechanism to produce gravitinos in thermal 
bath. If one adopts his mechanism, the constraints obtained below must become more stringent. However, 
it is not clear to us if his new mechanism is relevant. 
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photon scattering. But in the hmit of — 0, the total cross section for the photon-photon 
scattering is proportional to and this process loses its significance. Hence finally, photons 
are thermalized by (III) pair creation in nuclei, or (IV) Compton scattering off background 
electron. And (V) emitted high energy electrons and positrons lose their energy by the 
inverse Compton scattering off background photon. 

The Boltzmann equations for the photon and electron distribution function and /e are 
given by 



dt 
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(29) 



NEU 



where DP (PP, PC, CS, IC, and NEU) represents double photon pair creation (photon- 
photon scattering, pair creation in nuclei, Compton scattering, inverse Compton scattering, 
and the contribution from the v — v scatterings). In solving eqs. (|28|) and (^9]) , we assume 



fe{E) 



dt 



E 



UE) 



NEU 



dt 



+ 



UE) 



dt 



(30) 



which respects the energy conservation although the actual energy distribution may be differ- 
ent. This assumption is adequate for our purpose since the photon (and electron) spectrum 
is determined almost only by the total amount of the energy injection |T6|. Furthermore we 
treat muons (and anti-muons) as electrons with the same energy and neglect the contribution 
from tau leptons whose creation rate is much smaller than the other charged leptons. Full 



details for other terms are shown in ref. 16 



4 BBN and Photo-Dissociation of Light Elements 

The presence of gravitino might destroy the success of BBN. The gravitino may have 
three effects on BBN. First the energy density of the gravitino at T ^ IMeV speeds up the 
cosmic expansion and leads to increase the n/p ratio at the neutron decoupling time and 
hence ^He abundance also increases. Second, the entropy production due to the decay of the 



6 



gravitino reduces the baryon-to-photon ratio and results in too baryon-poor universe. Third, 
the high energy photons induced by the decay of the gravitino destroy the hght elements 
(D, ^He, ^He). Among three effects, photo-dissociation of the light elements by high energy 
photon is the most important for the relatively light gravitino {m^j2 < ITeV). The entropy 
production is not important for the present case since gravitinos decay into invisible particles 
(i.e. neutrinos). 

When the energy of the high energy photon is relatively low, i.e. 2MeV < £7 < 20MeV, 
the D, T and ^He are destroyed and their abundances decrease. On the other hand, if the 
photons have high energy enough to destroy '^He, it seems that such high energy photons only 
decrease the abundance of all light elements. However since D, T and ^He are produced by 
the photo-dissociation of ^He whose abundance is much higher than the other elements, their 
abundances can increase or decrease depending on the number density of the high energy 
photon. When the number density of the high energy photons with energy greater than ~ 
20MeV is extremely high, all light elements are destroyed. But as the photon density becomes 
lower, there is some range of the high energy photon density at which the overproduction 
of D, T and '^He becomes significant. And if the density is sufficiently low, the high energy 
photon does not affect the BBN at all. 

From various observations, the primordial abundances of light elements are estimated [ pUj] 

as 

0.22<rp=(^) <0.24, (31) 



PB 



p 



^"l >1.8xl0^^ (32) 
!^R±^) <1.0xlO-, (33) 

where p4He and pb are the mass densities of ^He and baryon. The abundances of light 
elements modified by the gravitino decay must satisfy the observational constraints above. 
In order to make precise predictions for the abundances of light elements, we have modified 
Kawano's computer code |^ to include the photo-dissociation processes. 

In the present calculation there are at least three free parameters, i.e. the mass of the 
gravitino ^3/2, reheating temperature Tr and the baryon-to-photon ratio rjB- However as 



shown in the previous paper |]I6[|, the baryon-to-photon ratio r]B is not important parameter 
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because the allowed value for tib is almost the same as that in the standard case (i.e. without 
gravitino). Therefore we fix r^^ = 3 x 10~^° in this letter. 

5 Results and Discussions 

The allowed regions that satisfy the observational constraints (pTD-(p^) is shown in the 
1^3/2 — Tr plane in Fig.|^. In Fig.^ one can see that for the gravitino with mass between ~ 
lOOGeV and ~ ITeV, the overproduction of D and ^He gives the most stringent constraint, 
while the upperbound on the reheating temperature is determined from the destruction of D 
for m3/2 — (1 — 3)TeV. Notice that D destruction is not considered in the previous work |[T3| . 
Furthermore the constrain from (D+^He) overproduction is more stringent. The reasons why 
we obtain the more stringent constrain are i) the gravitino abundance is (4-5) times higher 
than the one obtained in ref.|^ that the previous authors used, ii) the secondary neutrinos 
are taken into account in our calculation, and iii) the gravitino lifetime for t/'^ — > z/ + z/ 
is longer by a factor 2 than the lifetime for t/;^ — > 7 + 7. (In ref. it is presumed that 
r3/2(^M -^U + V)= r3/2(^^ ^7 + 7).) 

In Fig.|^, the overproduction of ^He due to the mass density of the gravitino at T ~ IMeV 
is also seen. The constraint from this effect is less stringent than those from the photo- 
dissociation of D and the present sneutrino mass density (discussed below). 

If we assume that the sneutrino is the lightest SUSY particle, the sneutrinos produced 
by the gravitino decay are stable and their present mass density is given by 

where h is the Hubble constant in units of lOOkm/sec/Mpc. Therefore the present limit on 



sneutrino mass 41GeV p2[ sets the upperbound on the reheating temperature; 

Tr < 6.6 X 10"/i^ GeV. (35) 

This constraint is also shown in Fig.^.0 

^In calculating eq.(^), we have ignored the effect of the pair annihilation of sneutrinos (i> + f + f)- 
For the pair annihilation process, the cross section is roughly estimated as av < aw/nT-f, (where aw is 
the coupling factor and v is the relative velocity), and the condition for sufficiently large pair annihilation 
rate {npav > H, with nc, being the sneutrino density) reduces to ni>/n^ > 10~^ (with Yj, being the yield 
variable for the sneutrino), which is less stringent than the constraint ( psf). However, if the sneutrino mass is 
small enough to hit the pole of the Z-boson propagator, the constraint (|35|) becomes weaker, and constraint 
from the ^He overproduction may become significant for large gravitino mass ^ 3TeV). 
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In summary, we have investigated the decay of the gravitino into a neutrino and a sneu- 
trino, in particular, the effect on BBN by high energy photons induced by the decay. We 
have found the constraint on reheating temperature and mass of the gravitino as 

Tr < 10^°"^^GeV for mg/s = lOOGeV - lOTeV, (36) 



which is more stringent than the one obtained in ref. []T3|, but compatible with the rough 
estimation of the reheating temperature (|l|) for a chaotic inflation. 
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Figure Captions 



Figure 1: Spectra of high energy electron-neutrino (sohd curve), muon-neutrino (dashed 
curve) and tau-neutrino (dotted curve) at T = leV for 777,3/2 = lOOGeV, ITeV and lOTeV. 
We have taken Tr = 10^°GeV. 



Figure 2: Constraints on the reheating temperature from BBN, the present mass density of 
the sneutrino. The region above the curves is excluded. The lower mass of the gravitino 
from 7773/2 > mp > 41GeV is also shown in long-short dashed hne. 
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